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its logarithmic partner, sourced by ordinary massless and by logarithmic non-normalisable 
gravitons, respectively. In addition, the logarithmic gravitons source two ordinary oper- 
ators, one with spin-one and one with spin-zero. The one-point function of the stress 
tensor vanishes for all Einstein solutions, but has a non-zero contribution from logarith- 
mic gravitons. The two-point functions of all operators match the expectations from a 
^ three-dimensional logarithmic conformal field theory. 
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1 Introduction 

During the last couple of years, a lot of effort has been put into the understanding of 
higher-curvature gravity theories in three space-time dimensions. In these studies, one 
centre of attention has been so-called "critical tunings" in the space of coupling constants. 
The corresponding theories exhibit several features making them potentially interesting as 
candidates for models of quantum gravity. 

The first instance of this criticality phenomenon was unravelled in cosmological topo- 
logically massive gravity (TMG) [1, 2]. The action of TMG consists of an Einstein-Hilbert 
term, a cosmological constant and a gravitational Chern-Simons term. In a seminal pa- 
per, Li, Song and Strominger [3] noted that there is a certain tuning of the Chern-Simons 
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coupling /i where several interesting phenomena simultaneously occur. First, one of the 
central charges of the boundary CFT vanishes and second, the linearised equations of 
motion degenerate. 

The fact that the central charge is zero means that, if the theory is unitary, it has to be 
chiral. Therefore, the partition function would be trivially holomorphic, fulfilling the hope 
[4] apparently not realised [5] for pure AdS gravity. This was taken as a strong indication 
that TMG at the critical point may have a consistent quantum mechanical formulation. 

In the same vein, the degeneration of the linearised equations of motion was taken 
as evidence that the massive graviton — a negative energy state for general tunings - 
is absent in the critical case. However, soon after the chiral gravity conjecture, it was 
realised that even at the critical point there is a physical bulk mode with negative 
bulk energy [6]. As in logarithmic CFT [7], the Hamiltonian was shown to have a Jordan 
block structure on the Hilbert space including ^ log and based upon this it was conjectured 
[6] that TMG at the critical point is holographically dual to a logarithmic CFT. 

The logarithmic mode has a different asymptotic behaviour than other modes, it di- 
verges more quickly. This is not a problem in itself — the variational principle is well- 
defined, and the required boundary conditions are consistent [8] . However, it opens up for 
the possibility to eliminate the mode by imposing certain boundary conditions 1 . It was 
subsequently shown [10], that the asymptotic behaviour cannot be avoided by constructing 
wave-packets with compact support because of a linearisation instability. 

There are therefore dual reasons to focus some attention on critical points in the study 
of higher-curvature gravity. On the one hand, they offer a potential mechanism eliminating 
problematic negative energy states. On the other hand — keeping the logarithmic modes 
— they offer holographic descriptions of logarithmic CFTs: an AdS/LCFT correspondence. 

As its discovery, the first technical pieces of evidence for this correspondence were 
provided in the context of TMG. Two-point functions were computed according to (suitably 
generalised) standard AdS/CFT methods [11], resulting in perfect agreement with LCFT 
correlators, and these results were generalised to three-point functions in a technical tour 
de force [12]. The one- loop partition function was also demonstrated to match LCFT 
expectations [13]. These developments have since been generalised to various other versions 
of higher-curvature gravity in three dimensions exhibiting critical points. (For an extensive 
list of references see, e.g., [14-16].) 

In particular, a certain combination of higher-curvature terms results in a theory that 
propagates unitary gravitons [17, 18]. This theory, which has acquired the increasingly 
misleading name "new massive gravity" (NMG), also allows for critical tunings leading to 
logarithmic behaviour [13, 19, 20]. 

Naturally, it is of interest to lift this entire discussion to higher dimensions. Such a 
venture was initiated by Lu and Pope [21] who formulated a higher-curvature gravity theory, 
being a four-dimensional analogue of critical NMG. This model was later generalised to 
arbitrary dimensions [22]. Subsequent works [23-25] have demonstrated and categorised 

1 We note in passing that a similar idea was advocated not long ago to eliminate the intrinsic ghosts of 
conformal gravity in four dimensions [9]. 
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the logarithmic excitations in the theory. However, the logarithmic structure is so far 
explored to a far lesser extent than in the aforementioned lower-dimensional literature. 
For instance, there are no results on any linearisation instability (potentially allowing for 
elimination of the logarithmic modes), there are no results on partition functions, and no 
correlators have been computed. 

In this work we take a first step toward a more complete higher-dimensional AdS /LCFT 
correspondence: we compute the one- and two-point functions of several operators in four- 
dimensional critical gravity. These operators are the stress tensor, its logarithmic partner 
and two other operators corresponding to the transverse vector and scalar pieces of the 
logarithmic graviton. All results are completely consistent with an LCFT3 dual. Our 
computations follow to a large extent the lead of [11] which also contains clear and com- 
prehensible explanations of the associated logic and technicalities. 

This paper is organised as follows. In section 2 we briefly review critical gravity in 
four dimensions and derive the first variation of the action. In section 3 we obtain the 
renormalised second variation of the action. This requires the construction and evaluation 
of several holographic counterterms. In section 4 we explicitly calculate the linearised 
modes and use them to derive the two-point correlators for critical gravity, and in section 
5 we conclude. Three appendices are devoted to technicalities that are not presented in 
the main text. 

2 First variation of the action 

To compute one- and two-point functions according to the AdS/CFT dictionary [26-28], 
one functionally differentiates the on-shell action with respect to the relevant sources. In 
this section we perform the first variation of the action and compute the corresponding 
boundary term. This object will be differentiated further in later sections. We also evaluate 
the first variation with respect to the boundary metric allowing to determine the one-point 
function of the stress tensor. The computation of one-point functions in general requires 
holographic renormalisation, and the method we follow was developed in [29, 30]. (For 
reviews, see [31, 32].) 

Our model is the one presented in [21] with bulk action 



Such higher-curvature actions generally propagate a massless spin-two graviton, a massive 
spin-two field and a massive scalar [33, 34]. The two latter excitations are ghosts. If the 
parameters a and j3 are tuned as a = —3/3, the spin-zero excitation is absent. If furthermore 
(3 = —1/2 A the theory becomes critical: the black holes have zero mass and entropy and 
the massive graviton becomes logarithmic [21, 23, 25]. 
Varying (2.1) with respect to the metric g^ u produces 
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where 

EOM^ = -Q*" - E" u , (2.3) 



(2.4) 



where the tensor A^ v is defined as 

A^ = (l + 2PR)g flu + 2aR flu . (2.5) 

The tensors appearing in the equations of motion are furthermore 

1 

2 



iiv — Rpv — ^Rg^u + A^jy , (2.6) 



>.\ 1 ,,\o - - • - 1 



E^u — ^{R^\R V — t-R a R\a9fiu) + 2(3R(R flu — -.Rg^v 



(2.7) 



2 

We shall from here on only consider the case a = —3/3. To reduce clutter we also fix A = —3. 
Using the (twice contracted) Bianchi identity it is easy to show that, for a = —3/3, the 
tensor E^ u is traceless. Thus, taking the trace of the equations of motion then establishes 

i? = 4A = -12. (2.8) 

In particular the Ricci scalar is constant on-shell, eliminating several of the terms in 
Eq. (2.7). Using again the Bianchi identities the equations of motion simplify to 

Qy.v + Y RXaR ^^u ~ 24/3(iV + 3g pu ) - 3f3(DR^ + 2R apXu R aX ) = . (2.9) 

Note that these are the correct equations of motion only in vacuum. Coupling the theory 
to matter would of course require keeping the full result (2.6)-(2.7). 
Let us now fix Gaussian normal coordinates, 

ds 2 = dp 2 + jij(p) dx i dx j . (2.10) 

From (2.2) it is then clear that the on-shell variation reads 



WbuiklEOM = i / d^xy/^yJf, (2.11) 



with JP defined in Eq. (2.4). Partially integrating the Christoffel symbols in J p allows us 
to write the variation as 

WbuiklEOM = A / d 3 x ^=7 J p = / d 3 x^f(- [A* + APP^SK^ 

J;>.\l -h)M v , 2 12) 



+ 



^V p (A ij + APP~f ij ) + V k A pk f j - VM« + A pp K ij 



This is the expression that we shall differentiate further to obtain the two-point correlators 
in later sections. In doing so, we shall use Poincare patch AdS4 as a background, to obtain 
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the correlators corresponding to a CFT3 on a flat background. This computation will also 
give us the one-point functions on a flat background. 

Let us first, however, compute the one-point function of the stress tensor in global 
AdS. The advantage of this is that we can then obtain the conserved Poincare charges of, 
e.g., a black hole. Our background metric is thus global AdS,i: 

ds 2 = -cosh 2 pdt 2 + dp 2 + smh 2 pdQ| , (2.13) 

and the matrix 7^ in Eq. (2.10) is assumed to have the expansion 

7« = 7gV P + 1% ~ tifpe-o + lfe-o + . . . (2.14) 
with the leading contributions fixed 

(-1/2 



(o) 

lij 



/-1/4 



V 




01/4 1 . 7 Jf 




-1/2 . (2.15) 



V 



-^sin 2 ' 



(3) (3) 

The term 7^ corresponds to the massless graviton and the term j5\- is forced to vanish 
by the equations of motion except at the critical point (3 = 1/6, where it captures the 
logarithmic mode. For other tunings, the massive graviton has a different power law fall- 
off. We ignore such terms for brevity, but still keep f3 arbitrary since it illuminates some 
of the results. It is important to keep in mind though, that our result for the stress tensor 
below is incomplete unless /3 = 1/6. 

Using these expansions it is straightforward to obtain an expansion for the variation 
in Eq. (2.11). This computation is detailed in appendix A and the result is 

WbuiklEOM = A / d 3 x V=7 ( [1 - 6/3] (K^djij - 2^5^) - ^ e -^g «5 7ij 

lK JdM \ 1 y ' 

(2.16) 

Here Kij is the extrinsic curvature 

Kij = \dplij ■ (2.17) 

(3) (3) 

The first term in this variation is divergent and contains variations of 7^ and , de- 
stroying a well-defined variational principle. These terms must in general be cancelled 
by holographic counterterms. Note however, that at the critical point /3 = 1/6 the first 
term vanishes. Thus, there is no need to add holographic counterterms for the critical 
case. This is in complete analogy with what happens in three dimensions for new massive 
gravity [20, 35]. 

Off the critical locus, counterterms are needed. It is shown in appendix A that the 
required terms are exactly the ones for pure Einstein gravity [36, 37] 2 , multiplied with the 
appropriate prefactor: 

hM = -^r I d 3 x v^T (4 - 2K + R[i\) . (2.18) 



2k 2 



dM 



2 For a less directly comparable, but earlier, computation of the counterterms for Einstein gravity, see 
[29, 30]. 
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The boundary stress tensor T 13 is denned as 

WrenlEOM = \f d 3 X T ij fry® , (2.19) 
1 JdM V 

or, shorter, 

T„ = -— J = ^ rGn|EOM , (2.20) 

where I ren = -^buik + ^a/Vf- At the critical point the stress tensor takes the simple form 

T crit — — /? (3) (2 2\) 

13 4k 2 ij ' 

while for generic values of j3 the result is 

^ = (1-6/3)^7^- (2-22) 

The latter result of course contains contributions from the boundary action in Eq. (2.18), 
and, for /3 = 0, is consistent with [37]. We note again that the result (2.22) only contains 

(3) 

the contribution from solutions captured by 7^ , e.g., massless gravitons and black holes. 
As shown in appendix A, the asymptotic equations of motion imply 

= 4)^ = ^ (3) ^ = > ( 2 - 23 ) 

and 



V 



k fl (3) 



(0)C ; =0, (l-6/3jV^=0, (2.24) 

where the covariant derivative V(o) is taken with respect to the boundary metric 7^ . These 

equations imply that the boundary stress tensor — for general /3 7^ 1/6 at least the part we 

(3) 

computed — is traceless and conserved. Note that 7^ is not necessarily transverse for the 
critical case. The non-transverse components make up part of the logarithmic graviton. 

Our result (2.21) shows that any solution having a vanishing p.-' also has a vanishing 
stress tensor. In particular this means that any solution to Einstein gravity has vanishing 
mass and angular momentum, confirming the result of [21] for the mass. The situation is 
completely parallel to that in NMG 



3 Second variation of the action 

To be able to compute the two-point correlators we first need the second variation of 
the action, or, equivalently, the first variation of the one-point functions. This entails 
computing the first variation of the coefficients multiplying 5"fij and 5Kij in (2.12). To 
technically simplify this computation we shall perform it perturbatively around Poincare 
patch AdS as opposed to the global metric. This means that the boundary is a plane, and 
that we obtain planar CFT correlation functions. 
Thus, we consider a metric of the form 



ds 2 = -y- + dx 1 dx 3 = — ^ h hij dx 1 dx 1 , (3.1) 
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where T)%j is the flat Minkowski metric on the boundary M 2,1 . The perturbation hij is 
assumed to have a Fefferman-Graham expansion of the form 

^ = bf^f + ^ + $ log. + $ + bf y log, + 4? y + • • • (3.2) 



To compute the correlator of the logarithmic mode we must include the leading logarithmic 

h(°) 



term that breaks the asymptotic AdS property of the metric. By the equations of 



(2) 

motion, this also requires the inclusion of the term logy. Furthermore, the equations 

of motion require b\^ be traceless. 

In AdS/CFT language, the expansion coefficients hfj and b^ represent sources for 
different operators. The former is a source for the stress-energy tensor (Tij), whereas the 
latter contains sources for its logarithmic partner, but also for other operators. Collec- 
tively, we shall denote these operators by tij. The tracelessness of b^ carries over to the 
operator(s) ty and imposes the constraint that tij be traceless too. Correlators of these 
operators are then given by functionally differentiating the on-shell action with respect to 
hf^bf. 

Our normalisation of these operators will be defined by 

( Tij ) = 2 5I ^p^ , (Tij ...) = 2/ J. (...) , (3.3a) 

(%)=2^jpi J ( tij ...)=-2i *(...), (3.3b) 

where the ellipsis denotes any operator. The factor of —i in the two rightmost expressions 
comes about because the generating function is actually Z-^ren (for Lorentzian signature), 
as explained in appendix B of [11]. 

3.1 Second variation of the bulk action 

The computation of the second variation of the bulk action is lengthy but straightforward. 
The details and the conventions are presented in appendix B. The final result is 



zon (0) oo ij 2 (0) *3 



8M (3.4) 



where all indices are raised by rf 3 . We also have already symmetrised 3 in the two variations, 
since only the symmetrised version is needed in computing the correlators. 

We see that setting b[j and o?) to zero recovers our old result (2.21) for the stress 
tensor at the critical point. Furthermore, because of the many divergent terms it is clear 
that the computation of two-point correlators requires holographic renormalisation, a topic 
to which we now turn. 



3 This means, e.g., setting S^^f +S 2 b^ 0) S 1 b\f = 26b^ 0) 6b\f and setting Sib^S^f -Ssbf^b^ = 
and so on. 
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S = I d 4 x^ 



M 



(2A/3 + 1)R - (2AB + 1) 2A + T F» U G^ + ——(F^F^ - F ) + 3F 



3.1.1 Auxiliary field formalism 

To determine the correct counterterms it is useful to consider first the formulation of a 
well-defined boundary value problem with the cut-off radius not taken to infinity. Since 
we are dealing with a higher-curvature theory, it is not enough to fix the metric at the 
boundary, but also some derivatives must be held fixed. A convenient and enlightening 
way to set up the problem, and to choose a combination of derivatives to hold fixed, is 
through the auxiliary field formalism developed in [17, 18, 35]. We shall use a minutely 
tweaked version, the tweaking being the non-linear version of the field redefinition used, 
e.g., in [19]. The field redefinition is such that the auxiliary field vanishes for the AdS 
vacuum. 

Running the risk of some redundancy, we present it in some small detail for general A 
and /3. We consider the bulk action 

-F^G + — 

(3.5) 

where F^ u is the auxiliary field, F = g^F^ and G^ u = R^ u — (R/2)g fMU . Varying with 
respect to F^, and using the trace of the resulting equation of motion, gives F^v in terms 
of the metric: 

BA 

FfMu = ^- ( — 6iV + R 9flu + 2A 9tll/ ) . (3.6) 

From this equation it is easy to see that F^ u = if and only if G^ v + Ag^ = 0, i.e., the 
auxiliary field vanishes exactly if the cosmological Einstein tensor vanishes. 

Substituting equation (3.6) into the action (3.5) yields after a little algebra the action 
(2.1) for a = -3/3: 

S = J d A x ^g [R-2A- spRm^ + PR 2 ] . (3.7) 

M 

Note that the action (3.5) is particularly simple for the critical tuning 2A/3 = —1. In fact, 
with A = —3 and B = 1/6 we have 



Scfl = Jd 4 xV=9 -F^G^ + \ (F^F^ - F 2 ) + 3F 

M 



(3.8) 



Apart from the term linear in F, this action looks precisely (up to a rescaling of F^ u ) like 
the auxiliary field formulation [18] of the pure higher-curvature theory studied by Deser 



The boundary value problem is now defined by requiring that the variation of g^ and 
Fpu both vanish at the boundary. Thus, the first variation of the on-shell action is allowed 
to contain boundary terms multiplying 5g^ u and 5F^ U , but no variation of the extrinsic 
curvature. Eliminating such terms requires adding a generalised Gibbons-Hawking term 
[35] which for the critical case (and in Gaussian normal coordinates) reads 

/ggh = ^ / d3x >/=7 F ij {Klij - K l3 ) . (3.9) 

dM 
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We have now set up a natural-looking boundary value problem for this theory: we keep the 
metric fixed at the boundary, and, since non-zero F^ v corresponds to non-Einstein modes, 
we keep the massive graviton fixed at the boundary. Note that this is nevertheless a choice 
— there are of course other possible boundary conditions. 

Fixing this choice now limits the number of additional allowed boundary terms. Most 
importantly, a boundary term may not change the boundary value problem. This limits 
us to use combinations of boundary intrinsic metric quantities and the auxiliary field Fij. 
In the next subsection we present a set of boundary terms that makes the second variation 
of the action finite. 

3.1.2 Renormalised second variation 

To regularise the action (3.4), one needs to find admissible counterterms that cancel all 
divergences. Since the computations themselves are not very illuminating, we defer them 
to appendix B. Let us here make some general comments. 

First, the generalised Gibbons-Hawking term (3.9) on its own is very complicated 
when expanded to second order. In particular, there is a non-zero contribution at the 
order ~ logy/y 3 . This is even more divergent than the terms present in the variation 
of the bulk action. Furthermore, the first variation of /qgh does not vanish. There is, 
however, a term that remedies these deficiencies. The term 




(3.10) 



dM 



where F = 7* J ' F^ can be used for this purpose. In fact, the combination 

I GG n-2I p (3.11) 

has vanishing first variation, and the second variation starts only at order ~ 1/y 3 - The 
full result for this term is presented in (B.24). This term alone, however, far from does the 
job. It only cancels one of the problematic terms in (3.4), and only at the price of adding 
another divergent term (at order ~ log y/y). 

To obtain a finite second variation one must add a number of terms involving the field 

(3) 

F^v, as well as the Ricci curvature R\- of the boundary metric jij. We have not been 
able to construct a finite action using only such terms however. Instead, to reach our goal 
requires adding a term of the form 

dM dM 

This term is Lorentz invariant at the boundary, but uses not only 7^ and Fij for its 
definition, but also the perturbation hij explicitly. (Or — equivalently of course — it 
uses the background metric 77^.) Although the counterterm (3.12) does not appear to be 
generally covariant we find that the boundary stress tensor is conserved. Thus the theory 
does not suffer from a diffeomorphism anomaly [40]. This suggests that the term (3.12) is 
actually equivalent to some covariant counterterm to second order in the perturbation. 
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The fully regularised action having a finite second variation reads 

Fen = ibulk + -^GGH ~ 2Ip 

- At f d 3 x J^i l-F^Fa - F ij df + —F ij V 2 F ij + —F ij (D 2 F) ii + hF . 
2k 2 J v ' L6 3 13 18 18 V 7 J J 

(3.13) 

Here, the differential operator D 2 is defined in (B.7) and computes, up to a factor, the lin- 
earised Ricci tensor around three-dimensional Minkowski space. Expanded in the Fefferman- 
Graham expansion (3.2) the second variation of I ren is 



(3.14) 



dM 

Already from this expression it is clear that all two-point correlators involving only Einstein 

(n) 

modes vanish: if all the 6y- are zero, so is the second variation. 
3.2 One-point functions 

The result in Eq. (3.14) also contains the result for one-point functions around a flat 
background. For our operator defined in (3.3a) we get 

whereas for the operator ty defined in (3.3b) we have 

3 ^(3) , 3 .( 3 ) 



<*«> = 2rf(5 + 2 9 * I (3 - IGi 
4 Two-point correlators 

The two-point correlators are given by the second variation of the action, or alternatively, 
by the functional derivative of the one-point functions (Ty) and (ty) with respect to the 
sources hffi and by . Differentiating the expressions in Eqs. (3.15) and (3.16) for the 
one-point functions according to Eqs. (3.3a) and (3.3b) we obtain 

(T ij (x)T kl (x')) = 0, (4.1) 
\ V K ) klV ,1 2 K 2 Slf 0) (x>) 2k 2 5h^(x) ' 1 ' 

{tij(x)t kl (x )) = -- 2 - jr _ + - -y— r . (4.3) 



^fl&j&CaO 2(56*' 



a:' 



Note that there are two ways to compute the (Ty ij^) correlator — either by differentiating 
(Tij) with respect to b$L, or by differentiating (tja) with respect to h l Ly In obtaining the 
expressions for the correlators we used the fact that 

OAl (0) 
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meaning that does not source any modes with logarithmic behaviour. We shall see 
explicitly that this is true when studying the modes below. 

As remarked earlier, correlators including only Einstein modes are identically zero. 
We immediately note that any different — log or otherwise — behaviour of the log-log 

(3) 

correlator with respect to the log-Einstein correlator stems from the g\- mode, as can be 
seen from Eq. (4.3). A logarithmic mode is only defined up to an arbitrary shift of the mode 
by Einstein modes, and this freedom can be used to eliminate the first term in Eq. (4.3). 
Thus, the non-trivial information in Eq. (4.3) comes from the last term: 

2 9i 5g {3) 

(tijtkl) = -^T ijtM )-^-^-. (4.5) 

Eliminating the first term would amount to redefining 

tij y t{j + g-^ij • (^ - 6) 

Thus, in order to compute the correlators, all we need to do is find the functional relations 
between the matrices b$ , h$ , b^ and g& . To achieve this we shall find all linearised 
modes in momentum space. This is the task of the next subsection. 

4.1 Modes 

(3) (3) 

We now aim to determine how the subleading Fefferman-Graham coefficients b\- and g\- 
functionally depend on b^ and nP. . This dependence comes about as a combination of 
the equations of motion and the boundary conditions in the interior of AdS. In the global 
case, the latter consists of requiring regularity, and in the present case of Poincare patch 
AdS, we require infalling boundary conditions at the Poincare horizon. 

In transverse gauge, the linearised equations of motion are rather simple. For the 
critical case they read 

(□ + 2)(D + 2)Vv = 0. (4.7) 

Here □ is the wave operator on AdS4. The solution space consists of Einstein modes tfi^ u 
and logarithmic modes ip 1 ^ satisfying 

(□ + 2)^ = 0, (4.8) 

and 

(□ + 2) 2 V^ g = 0, but (0 + 2)^^0, (4.9) 

respectively. Now, we need the modes expressed in Gaussian normal coordinates, which 
is not compatible with transverse traceless gauge in general. We shall anyhow proceed by 
solving (4.7) and then transform the solutions to Gaussian normal coordinates. 

To find the logarithmic, as well as the Einstein solutions, we solve the equation 

(□ + m 2 )?/v(m) = 0. (4.10) 
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Then the Einstein and logarithmic modes are obtained as 



and 



log _ d t 



III/ 



fJtV 



dm 



--V2 



(4.11) 



We work in three-dimensional momentum space, making the separation ansatz 



-ip-x 



(4.12) 



Using Lorentz invariance of the background, we can fix a certain Lorentzian covector p^. If 
Pi is timelike we choose p% = E 5j and if it is lightlike we choose pi = E Sj + E Sf . 



Timelike modes 



Let us start with the timelike case. Solving (4.10) with ip^ u (m) = e~ t/jfiv{pi',y), using 
the gauge condition V^"ifj^ u (m) = 0, is fairly straightforward. The most general solution 
has ten undetermined coefficients. 

An example of a solution where eight of the coefficients have been put to zero is 



e- lEt [C 1 UEy) + C 2 yu(Ey)} 



/o o o\ 

10 
0-10 
\0 0/ 



(4.13) 



where v = l/2(— 1 + iy/—17 + 4m 2 ), j v and y v are spherical Bessel functions and C\i are 
constants. (We shall present all solutions only in their final form.) 

The next step is to require infalling boundary conditions at the Poincare horizon at 
y = oo. To achieve this we note that 



[Ci j u (Ey) + C 2 y u (Ey)\ ~ C x cos(Ey 



---) + C 2 sin(Ey-- 



(4.14) 



from which it is clear that infalling boundary conditions correspond to C 2 = iC\. For 
the other solutions the reasoning is identical. In this way the number of undetermined 
components reduces to five. 

The next step is to construct six Einstein and five logarithmic modes by using (4.11). 
Four of the Einstein modes are pure gauge, whereas all the logarithmic modes are physical. 
The last step is to go to Gaussian normal coordinates. The mode (4.13) is already in 
this gauge. For those that are not, it is simple to construct the corresponding gauge 
transformation. 

Below we present the modes that are the result of these computations. 



Logarithmic modes 
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The full set of logarithmic modes is: 



logl 

ij 



iEtFi(Ey) 



log3 
ij 




log5 



-iEt 



/ 



'7 



F 2 (ify) 



V 



/-l o\ 
1 

V° V 



+ ^3(^2/) 




e- lht F 4 {Ey) 



y 



(4.15) 



(4.16) 



(4.17) 



The functions Fi(y) are defined by 
F 1 (Ey) = Ei(iEy)-e iE y , 



F 2 = Ei(iEy) - e iEy 



iEy 



F 3 (Ey) = (^X - 3)Ei(^y) + 3e 1 ^ 



F 4 (Ey) = e 



-iEy 



2iEy ^ + n(i + Ey) - 



(Ey - i)(ir + iEi(2iEy)) 



(4.18) 

(4.19) 
(4.20) 



It is clear that the first two modes correspond to a vector representation of the little group 



SO(2), that i^if 3, and ifl^ g4: correspond to a traceless tensor representation and that ifi^ e5 
corresponds to the scalar representation. It is therefore possible to immediately write down 
the modes for general timelike pi. 

1 2 

To this end we denote by e i ' an orthonormal basis in the orthogonal complement of 
Pi, and construct from them two traceless tensors M}- = eje] — e?e 2 and M}- = eje? + e]ef 
in the same space. Then, if \p\ = \/—p 2 , we have 



logl, 2 _ —ip-x Fi(\p\y) ( 1,2 



MA = e - ip 



If 



Pie 



1.2 



J 



Pj e i' 



2 ij 



-ip-x 



F2(\p\y)vij ~ F 3 (\p\y) 



PiPj 
\p\ 2 



(4.21) 
(4.22) 
(4.23) 



Einstein modes 

Regarding the Einstein modes, only the traceless tensor modes are not pure gauge. Thus, 
log-modes corresponding to scalar and vector Einstein excitations should be identified as 
Proca modes [25]. 

To find the pure gauge modes one simply makes an ansatz ifi^ = ^Gu£f) w hh = 
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e tEt £, u (y) and enforces Gaussian normal coordinates. The resulting modes read 



„— ip-x 

*f * = ( 4 - 24a ) 

,E3,4 = e -i(p- X -\ P \y) ^-i\p\y) M 12 (424b) 
e ~ip-x 

= — r- ( 2 ^- - > ( 4 - 24c ) 



^ 6 = ^Pipj ■ ( 4 - 24d ) 



e 



-ip-x 



Functional relations 

To derive how bf^ and 5^ depend on bff and /i-^ we must now construct modes that 
have either of these two expansion coefficients vanishing. The Einstein modes already have 
= 0, so they correspond to varying only . To get a mode that has vanishing a 
linear combination of the Einstein and logarithmic modes must be taken. To this end, an 
expansion around y = for the functions F^ is useful: 

F 1 (|p|y) = logy+[-l + 7 +y+log|p|]+^ + ^^ + ... (4.25) 
F 2 (\p\y) = logy + [ - - + 7 + - + log |p|] - + ... (4.26) 



(4.27) 



F 3 (\p\y) = - Slogy + [3 - 3 7 - — - 31og + logy 
+ ^[-3 + 27 + ^ + 21ogN] + M^ + ... 

F4(\p\y) =logy+ [-2 + 7-i7r + log2 + logH] + ^L^logy 

|/;|2//2 [ 7 -t7r + log2 + log|p|] -^j^logy (4.28) 



2 °" u 3 

i\p\ 3 y 3 



9 



■[-8 + 37 + log8 + 31og|p|] + ... 



Consequently we find the following combinations of logarithmic modes and Einstein modes 
with vanishing hff: 

(4.29) 
(4.30) 

(4.31) 



tit' 2 


= ^ + (1 " 


ITT 


- log |p|) ^g 1 ' 2 , 




= Vvf' 4 + (2 - 


7 + ITT — 


log 2 - log |p|) iff*' 4 


Ylj 




■ 27 — ITT 


-21og|p|) Vg 5 






- 7 "T 


-logN)^ 6 - 
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Let us start with the functional derivative 

x ( 3 ) 

This is the quantity of relevance for computing the two-point function of the energy- 
momentum tensor in Einstein gravity. Variation of g\- with respect to hjL. , keeping fel'L 

fixed, corresponds to analysing the relation between h 1 ^ and g\- for the Einstein modes. 
We let the index I run over 1, ... 6 and define the basis {e 1 } of symmetric matrices as 

{e 1 } = {fcjg} , (4.33) 

where the subscript hEl denotes the leading Fefferman-Graham coefficient of the corre- 
sponding Einstein mode tp EI from 4 Eq. (4.24). Similarly, by below, we shall mean 
the 0(y)-term in the Fefferman-Graham expansion of the corresponding Einstein mode. 
Explicitly we have 

ejj = Pie) + pje} , e| = pte) + pjef , (4.34) 

3_1122 4 _ 1 2 , 2 1 /a or\ 

e ij — e i e j e i e j ) e ij — e i e j ' e i e j > (4.60) 

e ij = ~ 2r lij , 4j = PiPj • ( 4 - 36 ) 

(3) 

The only modes having a non-zero are the traceless tensor modes (4.24b). Explicitly: 

E3,4 i\p\ 3 3,4 E^3,4 n (A nn\ 

9 m = -3~ e 4 ' = • (4.37) 

For a general Einstein mode hij = A-iwfj we therefore have 

h V = E M - 4f = M + M) • (4.38) 

To compute the desired functional relation, we only have left to compute how ^3 and A4 
depend on h$ , and to achieve this we need to invert the first relation in (4.38). 
Fortunately this is simply done. Noting that 

e^ej = 25 J K , K = 3, 4 and J = 1, ... 6 , (4.39) 

we have 

AK = \e%hf, # = 3,4. (4.40) 
Inserting this into the second relation in (4.38) we have 



4 



We actually take the Fefferman-Graham coefficient of the Fourier transform, i.e., we leave the expo- 



rt") .„,.] 

confusion, we shall write out the argument — x or p — explicitly. 



nential factor out. To avoid clutter, we do not put a tilde on h\- and g\J . Whenever there is risk of 
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and therefore 



(3) 



3 



on (o) u 

From the present form it is not clear that this expression is independent of the explicit 
choice of the polarization vectors e > 2 , but it is simple to show that this is the case. In fact, 
defining the matrix 

®ij(p) = VijP 2 ~ PiPj , ( 4 -43) 



and using that 



r lij = ^ + e}e} + e!e], {AAA) 



it is straightforward to show that 

S (3) 

5h^ = 6R (6 ^ 7 + @uQ > k ~ Ql ^ kl) ■ (445) 

This expression is identical to the Fourier transform of the two-point function of the stress 
tensor in a three-dimensional conformal field theory [41]. (See, e.g., Eqs. (105) and (108) 
in Ref. [42].) 

Let us now turn to the functional derivative 5b\^ /5b k ^, which follows from a very 
similar computation. Since we are differentiating with respect to b k L we are interested in 
the modes ip LI having vanishing h$ . We denote the corresponding Fefferman-Graham 
coefficients by bVL^ an< ^ so on ' 



Note first that we have 



1 k 3 c 3 



b m = e 1 , / = 1, • • • , 4 , b\l )l3 = —/- 3 e 6 = m - ^p iPj . (4.46) 

L p p 

For an arbitrary mode hij = Aj^f- 1 we furthermore have 5 



i=i 
and 



b V = E A A + ^ - p2P^) » ( 447 ) 



^ ) = -^(^4 + ^4) • (4.48) 
Using exactly the same construction as before, we therefore find 

5b {3) 5q {3) 

4 = -«& = -if (e "= e " + e>,e * - e « 8w) • (4 - 49) 



5 The source has only five degrees of freedom, because it is traceless by the equations of motion 
(B.15). For similar reasons gjjp and have only two degrees of freedom — they are transverse and 
traceless. 
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We have now only left to compute Sg^ /Sb^L. As noted before, the g^ corresponding to 
non-Einstein modes is not necessarily transverse (but always traceless), and these excita- 
tions correspond to several distinct operators in the CFT. Only the transverse traceless 
tensor part gives the logarithmic partner of the stress tensor. 

Using the York decomposition of a traceless tensor, we can split up the operator tif 

Uj = ViVj + VjVi + tj? + (ViVj - ^ V 2 ) S , (4.50) 

where V is a transverse covector {V l Vi = 0), is transverse traceless and S is a scalar 
operator. The logarithmic partner of the stress tensor is the operator tj^ . 

Correspondingly, we shall split up the (traceless) Fefferman-Graham components into 
three pieces: 

hij = ViVj + VjVi + hj^ + (ViVj - \riijV 2 ) s . (4.51) 



In Fourier space this translates to 



1 



with 



hij = -i (PiVj + VjVi) + hj^ - (piPj - -r]ijp 2 )s, (4.52) 

Pitf = and p^T = . (4.53) 

Thus, the transverse vector part of corresponds to the ejj and efj expansion coefficients, 
the transverse tensor part to e?- and efj, and the scalar part to the efj coefficient. 

Let us denote "transverse vector" by TV, "transverse tensor" by TT and "scalar" by 
S. Then, computing the derivatives 



tv' (56*' 



TT Sb^L 



(4.54) 



arbitrary mode h^ = Y^i Aiiph 1 to order y produces 



corresponds to letting only A\ t 2, ^3,4 and A$ be non-zero, respectively. Expanding an 

JLI 

\ 3 



+ M.( C - 2 log |p|) (A z e% + A4) (4.55) 



3 

'"'Vl ■ r~;< 



i\v\\ /I 5 , 3 „ 



18 ~~"V2~ ' p 2 
where 

C= 14/3- 2 7 + 7ri- 2 log 2 (4.56) 

is a numerical constant that might in the end be absorbed through redefining the loga- 
rithmic mode by adding suitable number times the corresponding Einstein mode. Note 
that the log \p\ term cannot be cancelled in this way. This is the term responsible for the 
logarithmic behaviour of the correlators in an LCFT. 
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To compute the desired variations, we now need to invert Eq. (4.47) and insert the re- 
sult in (4.55). Again this is straightforward, using that {ejj, . . . , efj, fcpfw} is an orthogonal 
basis, and 



ij 1,2 

e l,2 e ij 



-2N 2 , 



e ij e 3 ' 4 - 2 

e 3A ij ~ ■ 



jMij ,L5 _ r. 
b L5 b (0)ij ~ b • 



The result is 



with 



Sb to) 



X ( 3 ) 

tv <56^ 



X (3) 
TT <Sfc*j* 



(3) 



2z 

TV 9\p\ 



P(i&j)(lpk) 



Sg 



(3) 



Sg 



(3) 



TT ^(C-21og N )^ 



(3) 



PiPj 



P 



-Vij 



12 p 



P 



(4.57) 
(4.58) 

(4.59a) 
(4.59b) 
(4.59c) 



In (4.59a) the symmetrisations are taken over ij and M and are defined with the usual 
factor of 1/2. 

Spacelike modes 

The above construction of modes captures only the timelike case. There are no obstructions 
to a very similar analysis for spacelike modes — provided we make some obvious changes. 
For example we will have to choose different boundary conditions at the Poincare horizon; 
we demand the modes be regular in the bulk. Finally we end up with results akin to (4.45) 
and (4.59b): 



Sg 



(3) 



1 



[Qik&jl + &ilQjk — ©y©fci) , 



X (3) 
5 9ij 



TT 



C — log \p\ 
6|p| 



(@ik@jl + @il@jk — @y@fcz) 



(4.60) 
(4.61) 



Lightlike modes 

Similarly we can solve (4.10) with the ansatz ^^{m) = e lE ( t+x ^ flu (pi;y) and the gauge 
condition V^t/wC 7 ™) = 0- However, the modes that we find are all power series in y. 

Therefore imposing boundary conditions, such as non-singularity of the modes in the bulk, 

(3) (3) 

kills one half of all solutions. Moreover, we kill all possible and tyj . Thus, the 
integral over the lightlike momenta does not contribute to the correlators. In fact, when 
actually performing the Fourier transform in the next subsection we will temporarily pass 
to Euclidean signature and the complication of having lightlike momenta will not play a 
role. For completeness, we present all lightlike modes in appendix C. 
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4.2 The correlators 

We have already derived all results needed for computing the correlators in momentum 
space — in fact, up to numerical factors, the functional derivatives in the above subsection 
are the momentum space correlators. What remains is to translate these into a configura- 
tion space form. We shall perform this Fourier transform in Euclidean signature and then 
continue back to Lorentzian signature. 

To illustrate the procedure, a general (x) can be written as [c.f. (4.38)] 

h^(x) = J^sj d 3 pe-^ £ A 7 (p)e& , (4.62) 

where the factor l/(27r) 3 is purely conventional. An Einstein mode sourced by this h\®\x) 
would have an 0{y) contribution 

= (2^8 / d 3 pe-^(A 3 (p)ef 3 + Ai(p)eJ) . (4.63) 
Inverting the Fourier transform and solving for in (4.62) gives 

4m(p) = \4\a(p) f dV e^'h^(x') . (4.64) 

Thus, we have 

9§\*) = j£pj d **f dVe«— 0^(44 + et^hf^x') . (4.65) 

Comparing with Eq. (4.42), we recognise the Fourier space functional derivative in the 
integrand. The configuration space derivative reads 

Therefore, for the correlators we obtain 

(W^m=~^f^e'' J ^(p), (4.67) 

<t«M<«(0)) = -fcwtHto)) - pip^ / d V pl *j|j(!>). < 4 - 68 ) 

To compute these integrals we use the formula 

1 f p ir P x 1 T(- — a) 

~ C(a) / d 3 P -— — , C(a) = —j, -^—>- , (4.69) 



\ X \ 2a K ' J * |p|2(3/2-a) ' V / 4 a^3/2 r ( Q 

and its generalisation 



I rln|x| 2 C'(a) 1 



C(a) L |x| 2a C(a) |x 



2a 



e^lnhf 

^ U|2(3/2-a) ' ^-' U J 
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obtained by differentiation with respect to a. The formulas (4.69) and (4.70) suffice to 
calculate all configuration space correlators. 

As shown in [41] the two-point correlator of a spin-two operator Oij in three dimensions 
is given by 

{O i:j (x)O kl (Q)) = ^ (~ (hklji + Iuljk) ~ ImjVkl) , (4-71) 
hj = raj + , (4.72) 

where A is a numerical constant. We shall, however, find it more convenient to use the 
form advocated in [42] 

(O -(x)O w (0)> = AA ijM ^ , (4-73) 

where 

a^jw = ~ (6i fc e j7 + e«e ifc ) - -e^ re w , (4.74) 

®ij = Wj ~ ViP . (4.75) 

Note that 6jj defined in (4.75), is the Fourier transform of 6^ defined in (4.43), used to 
express the correlators in momentum space. Therefore, performing the Fourier transform, 
keeping A^jm, is quite trivial and we only need (4.69) and (4.70) for a = 1. For the 
correlator (4.67) we get 

(T tJ {x) t H (0)) = AA ljM ^ , (4-76) 

with 

The tij operator contains, as explained before, three different pieces. The transverse 
traceless part corresponds to the logarithmic partner of the stress tensor. It is clear from 
Eqs. (4.59) that the correlators of the other two operators — Vi and S from (4.50) — have 
the usual form for a spin-one and a spin-zero operator: 

(V i (x)V j (0))=A v Q i ~, (4.78) 
(S(x)S(0))=A S j^, (4.79) 



\x\ 2 



where 



1 7T 1 3-7T 

Av = j2^^ and As = ~j2^2^- (480) 

Note that the two-point correlator of the spin-zero operator is negative. This is another 
manifestation of the non-unitarity of the theory and means that the theory would be non- 
unitary even if the logarithmic partner of the stress tensor would be absent. 
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Evaluation of the transverse traceless part of (4.68) yields, via Eq. (4.59b), 

<^»^ T (0)> =AA ijM lQg |X|2 + g + 27 " 2/3 , (4.81) 

generalising (4.73). 

The ambiguity of the log-mode with respect to addition of Einstein modes is evident 
by the linearity of the operator Aj,-^. Thus, we can again use the freedom of redefining 
tj? to get rid of the factor C + 2 7 - 2/3 via, tf? -> ^ T - (C/2 + 7 - 1/3)% 

We have now found all two-point functions of four-dimensional critical gravity. The two 
correlators (4.76) and (4.81) together with the result (4.77) for the quantity A constitute 
the main quantitative results of this paper. 



5 Conclusions 

This work is part of an effort to increase our understanding of higher-curvature gravity in 
four dimensions. In particular, we study the critical tuning of the coupling constants in 
order to determine to what extent the AdS/LCFT duality, discovered in three dimensions, 
extends to the four-dimensional case. As a first step in this direction we computed the one- 
and two-point functions for the critical theory. 

There are four operators in the boundary theory, which are categorised by their sources. 
The stress tensor 2L- is sourced by fti? and its one-point function is transverse and traceless 
and given by 

The operators , Vi and S are sourced by the corresponding components of , and the 
one-point functions are 

+ !^f TT ) , (5-2) 

(5.3) 
(5.4) 

where gj^ ,TT , and are defined by the York decomposition of 

$ = V^ (3) + V;Vf > + ,f TT + (V,V, - 1 % V 2 )5(3) . (5.5) 

Note that the one-point function of the vector operator satisfies a Ward identity of the 
form V' l {Vi) = 0. 

The fact that the stress tensor vanishes for all Einstein solutions is consistent with the 
result that the mass and entropy of black holes in the theory vanish [21]. 



= 


3 

2^2 




{Vi) = 


9 

4^2 




(S) = 


9 

4k 2 


si", 
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The non-trivial two-point correlators all match the expectations for a logarithmic CFT, 
Tij and tf^ forming a rank-two logarithmic pair. Explicitly, 

{T ij (x)T kl (Q))=0, (5.6a) 
(T tJ (x)t^(0)) = A^|^p , (5.6b) 



/ tt / \ tt/„m 1 67r - log,(\x\ 2 m 2 ) 



where we parameterised the freedom to add a multiple of an Einstein mode to the log- 
mode by a fiducial mass scale m in (5.6c). The two-point functions of V{ and S have the 
form expected for ordinary spin-one and spin-zero operators (4.78)-(4.80), and all "mixed" 
correlators vanish. The relatively simple form of the log-log correlator (5.6c) suggests a 
natural generalisation to arbitrary dimensions: 

^ T (x)^ T (0))« A^j^g, (5.7) 



M) 1, A * A A , 1 



where 

To explicitly derive this result, and to compute the constant of proportionality, could be a 
worthwhile exercise. 

The result (5.6) demonstrates that a lot of the story from three dimensions is repeated 
in the present case. There are, however, also differences. In three dimensions, the loga- 
rithmic graviton corresponds to two degrees of freedom, each playing the role of one chiral 
component of the logarithmic partner of the stress tensor. In the present case, this part is 
played by the transverse and traceless tensor modes. But now there are three additional 
degrees of freedom, parameterised by Vi and S above. These were called Proca modes in 
[25], and their normalisable representatives do not have logarithmic fall-off at the conformal 

(3) 

boundary (b- = 0). The spin-zero operator gives rise to a negative two-point function, see 
(4.79) together with (4.80). This is another indication of the non-unitarity of the theory. 

The Proca modes have zero energy by (5.1). It would be very interesting to know 
whether these modes are subject to a linearisation instability akin to that found in [10]. If 
they do show a linearisation instability, then all logarithmic modes can possibly be trun- 
cated by imposing boundary conditions, leaving a theory propagating only massless (and 
energyless!) spin-two gravitons. Otherwise, such a truncation is impossible. We note that 
there are other consistency requirements for a truncation to work. The boundary condi- 
tions must be consistent and the higher-point functions between truncated and untruncated 
modes must vanish. See, e.g., the discussions in [12, 43]. 

Another interesting extension of the present work would be to extend the analysis to 
the higher-derivative gravity models recently presented in [44]. Degeneration of multiple 
massive modes with each other and/or the Einstein modes might lead to gravity duals for 
higher-rank LCFTs in arbitrary dimensions. Such an extension would then involve the 
calculation of two-point functions in higher-rank Jordan cells. 
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Furthermore, while critical gravities are non- unitary [45] (as are all duals of LCFTs), 
critical tunings of higher-derivative models that lead to odd-rank LCFTs, might open up 
the interesting possibility of unitary truncations [46]. First steps of such a venture in the 
gravity context are taken in [47]. 

Last, but not least, it would of course be very rewarding to find a suitable condensed 
matter application for higher-dimensional LCFT, and to construct a phenomenological 
holographic model of it, using the model explored in this paper as canvas. 
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A Global coordinates: Fefferman— Graham expansion and first variation 

A.l First variation of the action 

In the Gaussian normal coordinates of Eq. (2.10) we have 



The matrix 7^ is assumed to have the expansion (2.14) with the leading contributions 
fixed. 

In Gaussian normal coordinates and after partially integrating the variations of the 
Christoffel symbols, the boundary term coming from the first variation of the bulk action 
becomes 





(A.l) 



5/bulk|EOM 



1 

2^2 



JdM JdM 




(A.2) 




For the case of interest a = —3/3, the tensor E^ u of Eq. (2.3) is traceless, and the trace 
of the equations of motion thus implies R = 4 A = —12. Asymptotically this implies the 



tracelessness of f3^ and 7®: 



^ /3 (3) - 7?)^ = ^ 7 (3) - 7?)^ 



0. 



(A.3) 



Below the symbol « denotes equalities after these equations have been imposed. A quantity 
of interest is the Ricci tensor. We obtain the following expansions for its components: 



Rij = -37ij + 



-37 



(-3p + l)Tr /3 (3) + 3Tr 7 (3) 



-2P\ 



R 



ip 



-e~ p R iy = 0(pe- 3 P) , 
R pp = e- 2p R yy = -3 + l -e^ p ([3p - 4]Tr /3^ _ 3Tr 7 < 3 )) + 0{pe~ bp ) 
« -3 + 0{pe~ bp ) , 

R » -12 + C(e- 4 ^) . 
Other useful expressions are (indices are raised by 7^; 



lJ o) = [7 ( - ;n 



7 



A", 



7(o) e 7(2) + A> e P( 3 ) e 7(3) j 



'J 



i^=^- e - 2 "7i 2)J +pe 



3 „ _,., r .l 



K* = e-^7fo) 



2e" 4 "7^ + 



7( 2) + 2 Pe 



>^(3) 



%j + «7(3)j] + 



e" 5 ' \f?> + 57^] + • • • 



2 " ^(3) 1 " '(3)J 

^ = (e 3 P -eP)+... 

We use these expressions to obtain the following expansions 

A ij « (1 - 6/3) 7 * j + 9/3 e~ 5p /3g } + . . . 

A pp = A" « (1 - 6/3) + O(e' 4 , 
« -9/3/3^ e" 3 " + 0(e- 4 ") , 
V fe ^ fc « 0( e - 4 ') , 
V p APP « 0(e- 4 ^) , 
V P A« « -27/3/3g )e - 5 ^ + ... 
Putting these expressions together yields for the variation of the on-shell action 



^bulk|EOM = I d a x 



1 

2^2 



dM 



d d x v^7 [1 - 6f3]{K lJ 6~f i:j - 2f 3 5K, 



27/3 



13 J 



-5poij 



(A.4) 

(A.5) 
(A.6) 
(A.7) 

(A.8) 
(A.9) 
(A.10) 
(A.ll) 

(A.12) 

(A.13) 

(A.14) 
(A.15) 
(A.16) 
(A.17) 
(A.18) 



(A.19) 
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which is the expression quoted in the main text. 

For P = 1/6 we do not need any holographic counterterms to get a well-defined varia- 
tional principle. This is similar to the three-dimensional case, where we did not need any 
counterterms for the logarithmic point of NMG [20, 35]. The stress tensor can be directly 
read off from the above equation and is proportional to . 

If /3 1/6 (in which case )3%\ vanishes by virtue of the equations of motion) we 
need a holographic counterterm. It turns out that the same counterterms as for pure 
four- dimensional gravity [36, 37] multiplied with a numerical factor does the job. 

After adding this counterterm, the action reads 



2k 2 
1-6/3 



\ [ d 4 x yf^j \R - 2A + MR 2 - ZWR 
J M l 

/ d 3 x^/=7(4-2K + i?[ 7 ] 

JdM V 



(A.20) 

d 3 x7=7(4-2K + i?[ 7 ]) , 

dM 



where R[y] is the contracted Ricci tensor of the induced metric 7™ on the boundary. 
Expressed in the components of the full Ricci tensor this quantity is given by 

fl[ 7 ]« = Rji + d p K ij + KK tj + 2K il K\ . (A.21) 

Using the expansions (A.4)-(A.ll) it is now straightforward to obtain 

R-hf = -l^f~ Ap + ^lf~ 4p ■ (A.22) 
Armed with this equation we derive the variation of the boundary term: 

(A.23) 



-7(4 -2K + R[ 7 })j = V=7( - n^Kij + \2K^ - 7 « + 7 { if- 4p ]6 7ij 
+ 0{pe- p ). 



This allows us to determine the variation of the full action (A.20). Adding the contributions 
from Eqs. (A. 19) and (A.23) produces = 0) 

WrenlEOM = ^ [ d 3 * | (1 - 6/3) 7 g) ^ • (A.24) 



WJdM V 2 

This shows that the stress tensor in this case is proportional to 7 A . As noted in the main 
text, because we ignore other fall-off behaviours than present in (2.14), this result does not 
contain contributions from massive gravitons, but only from massless gravitons and black 
holes. 



A. 2 Asymptotic equations of motion 

Let us now turn to the asymptotic equations of motion. By taking the trace of the EOMs 

(3) (3) 

we already established tracelessness of 7^ and ■ To show that the stress tensors for 
the critical case (/3 = 1/6) and the non-critical case (/3 7^ 1/6) are conserved we need to 
take a look at the ip-components of the EOMs. For /3 = 1/6 we use the ansatz (2.14) and 
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plug it into the equations of motion. Using the simplification that follows from that R is 

(3) (3) 

constant and throwing away terms proportional to the traces of 7^ and p^- we obtain 

EOM, p = ~ V^gV 3 ' + 0{pe-^) . (A.25) 

Here V(o) is the covariant derivative with respect to 7 |? • This proves the conservation law 
for the critical case. 

In the non-critical case ^ 1/6 we do not have logarithmic modes and therefore we 
(3) 

have to omit the term in the expansion (2.14). Again the ip-components of the EOMs 
yield 

EOM ip = ^(1 - 6/3) V^gV 3 " + 0(pe-^) , (A.26) 

up to trace terms. Thus, the stress tensors (2.21) and (2.22) for the action (2.1) are traceless 
and conserved for any value of /3. 



B Poincare coordinates: Fefferman— Graham expansion and second vari- 
ation 

We expand around Poincare patch AdS: 

d S 2 ''f • 7ii dr< d*> = ' "t 1 ''' 11 " + /',, «!.'•' <I.H , (B.l) 

yZ yZ 

where rjij is the flat 3d Minkowski metric, and h{j has the expansion 

, JO) logy (0) 1 (2), (2) . ,(3) , . (3) s 

hi i = b V~y2- + % 75 + b \j lo %y + % + y logy + % y + ■ • • ( B - 2 ) 

near y = 0. Our goal is to compute the variation of the terms multiplying the variations 
<5 7 jj and <5-fQj in (A. 2). To give them names, let us define 

MbuiklEOM = 7T2 f d3x J p = f d 3 x (tfjSKij + Rp-fij) . (B.3) 

Zk JdM JdM V 7 

We have computed and i? 2 J only on-shell in (A. 2), but this is all we need. In fact, 
to compute the correlators, we put in variations that satisfy the linearised equations of 
motion. Therefore, it is enough to vary the on-shell quantities, and we may, along the way, 
use on-shell relations between the Fefferman-Graham coefficients. 

For the purpose of expanding R£ and R% , and to expand the equations of motion, 
we need the Ricci tensor to the first order in hij. It is most simply obtained using some 
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computer algebra package (we used GRTensorll [48]) and it reads 

,2 

fnVl,; + — 4- Hop- 7/ - 1 , ,, 



R 



Kl 



o !/ /n 2 n ,3 (0) /n\J2) . (2) 37/ ( 3 ) 

3 7ij " y P 7)ij + ^2^/ + (logy - 1/2)6}/ + - —b\/ 



+ 



1 



Tr 6 (0) + (logy + l/2)Tr 6 (2) + Tr g 



(2) 



+ | (1 + 3 logy) Tr + ^Tr g ( 3 )] % - + 0(y 2 logy) + 0(h 2 ) , 
Riy =Y ( V ^(0), " V * Tr 6(0) ) + ( lo §y + V2) 2/ (V*6f 2)i - V,Tr 



2y V ■ ^ 
+ (Vfc5( 2)i " ViTr y( 2 )) , 



i? 



1 



yy 



_ 3 + Tr 6(°) - y 2 Tr b& - ^-(3 logy + 4)Tr 6^ - ^Tr :1 
y V 2 2 



(B.4) 



(B.5) 



(B.6) 



In these equations, and always when it comes to the quantities g^' and the trace is 

taken with respect to the metric rjij, i.e., Tr b^ = jf^b^ . Also, the indices of g^ and 

b^j are always raised by rfi and all covariant derivatives are with respect to rjij, i.e., are 
just ordinary derivatives. We also defined the differential operator D 2 as 

(2>V B) )« = ViVj-CTr 5 (n) ) + V 2 ^ - (V,V fc yf n)j + V,V fe yf n)l ) . (B.7) 

Note that the operator D 2 computes, up to a factor of —1/2, the linearised Ricci scalar 

(n) 

corresponding to a metric perturbation g\- around a flat background. 

B.l Linearised equations of motion 

Let us start by expanding the linearised version of EOM^ as 

oo 

EOM«= [(EOM^)^y n logy + (EOM n )^y"] 

n=-l 

and the linearised Ricci scalar as 



(Bi 



R^ = Y,[Rn.y n ^gy + Rny n } 



n=0 



(B.9) 



Evaluation of these quantities, which all must vanish, yields the linearised equations of 
motion. From the constancy of the Ricci scalar we obtain the following relations (i?g , R^ 
and R\ are all manifestly zero) 

(B.10) 
(B.ll) 
(B.12) 
(B.13) 
(B.14) 



R{) 


= - 4Tr 6 (0) 


= 0, 




Ri 


=V 2 Tr 6(°) - 


" ViV,-6j J 0) - 


- 4Tr b {2) = , 


R2 


=V 2 Tr - 


-Wv>f Q) 


- 4Tr y (2) = 




= - 3Tr 6 (3) 


= 0, 




R 3 


= - 3Tr y (3) 


- 2Tr b {3) 


= o, 
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from which we deduce 



TV 6 (0) = Tr # (3) = Tr 6 (3) = , 

ViVi&jJj+^&f 2 ) =0, 

V 2 Tr - ViV^j J 0) - 4Tr </ 2 ) = . 



(B.15) 
(B.16) 
(B.17) 



Using these equations to simplify the expressions coming from (EOMq )ij and (EOMo)ij 
yields 



( J DV°)) y =26g ) + 2Tr6( 2 ) % . 



1 -2^(0) o.(2) . „ (2) 



-V z 6 
2 « 



3 ¥ + w + 



2Tr 5 (2) + Tr 6 (2) 



(B.18) 
(B.19) 



B.2 Second variation of the action 

Using the equations of motion derived in the last subsection, we now want to compute 
the action to second order in the perturbation and expand it in y. Also this is best done 
using computer algebra. When we present the result we shall not keep track of the two 
variations separately since, for computing correlators, all we need are the symmetrised 
variations. More explicitly, if, say 



8^S[8ijij,5 2 ^ij} 



1 

2/? 



(B.20) 



dM 



we shall only write 



6®S 



2k 



dM 



(B.21) 



In this appendix, we shall often even leave the integral sign, the factors of 1/2k 2 and \/— 7, 
and the 5 out and write 



(B.22) 



when actually meaning (B.20), hoping that this does not lead to confusion. 

Using these conventions, and using the equations of motion, the second variation of 
the bulk action reads 



6®L 



bulklEOM 
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B.3 Boundary terms 

We now apply the same computational and notational framework to various counterterms. 
The terms in question are {F^ u is the auxiliary field and F = ^Fij) 

1 



IgGH 

Iff 
Ifvf 

hp 



i 
1 
1 
1 

2k 2 



dM 



d 3 x^F ij F i:j , 



dM 



d 3 x^/=iF ij V 2 Fij , 



It 



Ifr 



l FDF 



2k 2 



(Tx^--fF , 



dM 



2k 2 



dM 



2k 2 



dM 



dM 



d?xy/=^F i ^{D 2 F) i j, 
(B.23) 



dM 



The second variations of these quantities are 6 



+ 



(0) ij 



- 3 %*f " « h W + 6Tr b(2)Tt h(0) + 66 (o)4? 



= -3&M>^ + [-3J& &<? + 3Tr ^Tr fc<°> 



(0)"ii 



(0)"ij 



5^Ifvf 



-27&gj&g> - 18/^if + 18Tr ft^Tr /»«» + 18&g )fl g> 



^ 2 )/^ = 2Tr ft^Tr Z^-. 

y 

Combining these results and defining 

1 15 
-ftot = -fbuik + ^ggh — 21 p — qIff + Ifr — Tg^ FVF — 18^ FDF ~ ^~ hF ' 

it is straightforward to show that all divergent terms cancel. The final result is 



8&L 



tot|EOM 



1 

2^2 



d 3 X\/-7 



2 (0) ij ^ 4V"(0) y ij 



(0) ij 



(B.24) 

(B.25) 
(B.26) 

(B.27) 

(B.28) 

(B.29) 

(B.30) 
(B.31) 



9Af 



6 The generalised Gibbons-Hawking term 7ggh actually has a non-vanishing first variation, but a linear 
combination of 7ggh and Ip has vanishing first variation. The variation of this combination is displayed in 
(B.24). 
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C Lightlike Modes 



To obtain the lightlike modes we solve (4.10) with the ansatz ( m ) = e * E VwCPtll/) 
and the gauge condition V^tp^lm) = 0. Going to Gaussian normal coordinates and 
choosing combinations of log and Einstein modes such that all hf^ are zero we find the 
following log-modes: 

, /o l\ 

aogl _JE(t+x 1 ) c3 logy + c4y J logy 



(C.l) 



,\ og 2 _*E(t+x 1 ) c9 log y + clOy 3 logy 

Y/J,l/ e 



1 
V i 10/ 

/no\ 

110. (C.2) 

^000/ 

These are vector modes orthogonal to our chosen pj. Three of the other log- modes turn 
out to be Proca modes, their respective Einstein modes are pure gauge modes. Note that 
there are no explicit logarithms, but they are identified as log- modes by the equations of 
motion: 

A) l\ 

1 , (C.3) 

\iio) 



^Procal = ( 



^+* 1 ) C 2y(l + ^) 



Proca2 



JEit+x 1 ) 



c6y 



Finally, we also have 




(C.4) 



(C.5) 



lo g 3 _ JEit+x 1 ) 

V* y2 



log4 



JEit+x 1 ) 




(3 + 2 logy) 



(3 + 2 logy) 




(C.6) 



(C.7) 



and 



log5 

flU 



JEit+x 1 ) 



c5 r 



logy 



,2 E 2 



A o o\ 

10 

y0 0y 



+ y -—(2 + logy) 




y A E A 

I | +^-(5-41ogy) 




1 1 
\0 0, 



(C.8) 



(C.9) 
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The corresponding 'simple' Einstein modes read 



/ Einstl _ JEjt+x 1 ) °3 + I/ 3 c4 




110 

Einst2 =c ^q+x 1 ) c9 + y J clO [ 110 



,0 0, 



The Einstein modes corresponding to the log-modes ^ log3 , V' log4 an d V' Iog5 are 

/o 1 N 



Einst3 _ iEit+x 1 )^ 



, Einst4 _ iEit+x 1 )^ 
Yfiu — c o 



,/,Einst5 _ AEit+x 1 ) ^ 
"^V ~ C y 2 
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] 




U10J 





10 1+ y 2 £ 2 I 1 1 



2 F 2 




4p4 



y 4 £ 



v oooy y o o -iy 

The gauge modes should be constructed from 

(l 1 




Egaugel = JEit+x 1 ) ^_ 



ge3 _ JEit+x 1 ) ^_ 

'\W y2 



^Egaug 



,Egauge4 iBf't+x 1 ) jL 

e y2 




^Egauge2 = ^(t+x 1 ) ^_ 




- 1 o t f I 1 1 o 

,0 0, 



(CIO) 



(C.ll) 



(C.12) 



(C.13) 



(C.14) 



(C.15) 



(C.16) 



(C.17) 



As mentioned in the main text all lightlike modes are power series in y (some multiplied 

(3 s ) (3) 

by logarithms). Requiring non-singularity of the modes at y — > oo kills all the g\j and b\- 
parts. Therefore, the lightlike modes do not contribute to the correlators. 
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